Magnetic excitations in dipolar pyrochlore antiferromagnet Gc^S^O^ 
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The spin dynamics in the geometrically frustrated pyrochlore antiferromagnet Gd2Sn207 is stud- 
ied by means of the electron spin resonance. In the ordered phase (TV = 1 K), we have detected three 
gapped resonance modes. Their values agree well with the developed spin-wave theory which takes 
into account the Heisenberg nearest-neighbor exchange, the single-ion anisotropy and the long-range 
dipolar interactions. The theory also predicts a fourth lowest-frequency gap, which lies beyond the 
experimental range of frequencies, but determines the exponential decrease of the specific heat at 
low temperature. 

PACS numbers: 75.30.Sg, 75.50.Ee, 75.30.Kz. 



I. INTRODUCTION 

The geometrically frustrated Heisenberg antiferromag- 
net on a pyrochlore lattice has an extensive degeneracy of 
the ground stated As a result, weak additional interac- 
tions play a prominent role leading to a multitude of mag- 
netic phases and phenomena. A typical example is pro- 
vided by Gd2Sn2C>7 and Gd2Ti20j, two pyrochlore anti- 
ferromagnets with close values of the nearest-neighbor ex- 
change constants, which, nevertheless, have different or- 
dered magnetic structures. While the antiferromagnetic 
structure in Gd2Ti 2 C>7 is described by the ordering wave- 
vector q = (1/2,1/2,1/2) and has not yet been com- 
pletely resolved^ the second material Gd2Sn207 has a 
much simpler four-sublattice spin-structure with q = 0. 3 

An exponential decrease of the low temperature spe- 
cific heat of Gd2Sn207 was recently found by the calori- 
metric measurements^ This result unambiguously points 
at a gapped excitation spectrum in the ordered magnetic 
phase. It is somewhat unexpected in view of the re- 
sults of analogous measurements for Gd2Ti2C>7j 5 which 
yield a power-law behavior C oc T 2 down to temper- 
atures of the order of 100 mK. Presence of non-frozen 
magnetic degrees of freedom at very low temperature 
was also related to the persistent spin dynamics observed 
via the muon spin relaxation (/xSR) in Gd2Ti207^^ and 
via Mossbauer— and /xSR measurements^ in Gd2Sn2C>7. 
The observed difference in the low-temperature depen- 
dence of the specific heat in the two pyrochlore antiferro- 
magnets must be further corroborated by a direct investi- 
gation of the excitation spectrum. Since inelastic neutron 
scattering measurements are hindered in Gd-compounds, 
the most convenient experimental technique left is the 



electron spin resonance (ESR). An absorption of pho- 
tons with a typical wave-length A ~ 0.1 -j- 10 cm provides 
a high- resolution probe of low-energy excitations. For 
the most probable single photon-magnon process only 
magnons with small wave- vectors k ~ 1/ A —> can be 
excited. 

Our previous measurements of the ESR spectra in 
Gd2Sn207 were performed at temperatures between 10 K 
(the Curie- Weiss constant is |#cwH — 9 K) and the or- 
dering transition at TV = 1 K. 1 - This temperature range 
is characterized by strong spin correlations in the ab- 
sence of an order parameter, which is commonly called a 
cooperative paramagnet. An unusual transformation of 
an exchanged narrowed paramagnetic line into a single 
gapped resonance mode with linear field dependence was 
observed in this regime. 

The present work investigates further transformations 
of the resonance spectrum at temperatures below the 
magnetic ordering transition. Three gapped modes are 
observed, two of which are degenerate in zero magnetic 
field. The gap values are reproduced by the spin-wave 
calculations using the known values of the exchange con- 
stant, the dipolar and the single-ion anisotropy energies. 
The fourth gap predicted by the theory remains unob- 
served in our experimental range, but agrees well with 
the specific heat data. 



II. SAMPLES 

Powder samples of Gd2Sn2 07 were prepared by the 
method described in a previous publication^ For sam- 
ple characterization we have measured the specific heat 
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in a 3 He- 4 He dilution fridge in the temperature range 
0.1-2 K. The data are obtained using a quasi-adiabatic 
technique with continuous heating of the sample.— In 
order to improve the thermal contact between the insu- 
lating powder sample and the sample holder, 1.5 mg of 
the powder was mixed with 1 mg of Apiezon N grease, 
wrapped into a silver foil and pressed. The quality of the 
thermal contact has been verified by comparing the data 
obtained in experiments with different heating rates. The 
range of fully reliable results corresponds to temperatures 
above 0.15 K. 

Our experimental results for the specific heat are pre- 
sented in Fig. [TJ Two subsequent temperature scans 
shown by open and closed circles (slower and faster scans 
respectively) are in satisfactory agreement with each 
other. A sharp peak is observed at the ordering transition 
with Tat = 1.01 K. The previously published data^i ob- 
tained with the same measurement technique are shown 
by squares and lie 5-10% lower in the whole range of 
temperatures. The measurements from Ref. (triangles) 
demonstrate a perfect agreement with our data above 
0.4 K, where both curves can be empirically fitted by a 
C oc T 2 dependence. Both sets of data also agree, but 
less perfectly, below 0.4 K and show a drastic decrease on 
cooling. According to our measurements, the most sig- 
nificant reduction of the specific heat (by a factor of 20) 
occurs in the temperature interval 0.3-0.15 K. Decreasing 
temperature further down to 0.1 K leads to an additional 
specific heat reduction by a factor of 5 to the level, which 
does not exceed the experimental accuracy. The specific 
heat decreases faster than any reasonable power law T n , 
which suggests an exponential temperature dependence 
in agreement with the previous study^ 



III. MAGNETIC RESONANCE 

Magnetic resonance measurements have been carried 
out in a transmission type spectrometer with a cylinder 
cavity designed for frequencies above 25 GHz, which was 
equipped with a 3 He cryostat with a minimum working 
temperature of 0.4 K. The magnetic field up to 100 kOe 
is generated by a cryomagnet. The absorption spectra 
are recorded on forward and backward field sweeps. 

To start with, we have traced the evolution of the res- 
onance absorption lines on cooling the sample from the 
strongly correlated disordered state through the ordering 
transition at T/v = 1.0 K down to the lowest experimen- 
tal temperature (0.45 K), at which the system is fully or- 
dered (the magnetic transition is first order—) . Two sets 
of measurements with different polarizations of the mi- 
crowave field with respect to the external magnetic field 
were performed. The left panel of Fig. [2] shows the reso- 
nance spectra of a sample glued onto the bottom of the 
cavity where the microwave field has only a component 
perpendicular to the external field h mw _L H. The sin- 
gle resonance line observed at all frequencies for T > 
changes its shape and shifts to larger fields when going 
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FIG. 1: (Color online.) Thermal variation of the specific heat 
divided by temperature in Gd2Sn207 measured with differ- 
ent heating rates (• and O are slower and faster scans re- 
spectively); □ and A are previous data from Refs. HTll; the 
dashed line is a linear approximation to the high temperature 
part of the data, the fit by a solid line is described in the text, 
the arrow marks the low temperature limit of data reliability. 
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FIG. 2: Evolution of the magnetic resonance absorption spec- 
tra in a powder sample of Gd2Sn207 on cooling from above 
the ordering transition to T = 0.45 K, recorded in a perpen- 
dicular polarization of the microwave field with respect to the 
external field h mw -L H (left panel) and in a tilted polariza- 
tion (right panel); the recorded lines are shifted upwards for 
clarity. 



through the transition. This shift results from an un- 
usual linear field dependence of the resonance gap in the 
cooperative paramagnetic state, as discussed in Ref. ITol . 

The absorption spectra are significantly modified when 
the sample is placed into a microwave field with a com- 
ponent along the external field. Additional lines develop 
in the spectrum below 1 K, one of them having a much 
larger intensity than the others. The properties of all 
these resonance modes were studied in detail at the low- 
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FIG. 3: (Color online.) The magnetic resonance absorp- 
tion spectra in a powder sample of GcbS^O? recorded at 
T — 0.45 K for various frequencies with h mw _L H (upper 
panel) and in a tilted direction of h mw (lower panel); dashed 
lines are guide-to-eye for tracing the field evolution of different 
resonance lines labeled by letters A-E; the inset in the lower 
panel expands the in-frame part of the absorption record at 
v = 62.0 GHz including spectral lines D and E. The transition 
into a spin-polarized phase at H 3 — 53 kOe is marked by a 
vertical line. 



field H s = 53.0 ± 0.5 kOe. A careful study of the spec- 
trum above the critical field again reveals the existence of 
three weak components (line D and double split line E) 
increasing in field, with the line E doublet exhibiting an 
exact linear field dependence. 

Our ESR measurements directly detect three reso- 
nance modes in zero magnetic field. Two of them (lines A 
and B) are exactly degenerate with a gap A2,3 = 1.61 K 
and the third mode has a larger gap of A4 = 4.1 K. In 
addition, the previous specific heat measurements and 
theoretical calculations suggest the presence of a fourth 
lowest branch at energy Ai ~ 1.2 K4*i^ The correspond- 
ing frequency v\ ~ 25 GHz lies near the lower boundary 
of our experimentally accessible range. Hence, it is natu- 
ral that this mode remains unobserved in the ESR experi- 
ment provided that it decreases with increasing magnetic 
field. A crude estimate from the low-T fit of our specific 
heat data C(T) cx T" 1 / 2 e~ A/T (the corresponding fit is 
shown by the solid line on Fig.[T]) also gives Ai ~ 1.0 K. 
One should note here that more elaborate fits at tempera- 
tures below 0.15 K (including, e.g., nuclear contributions) 
would rely on an overestimation of the experimental ac- 
curacy due to the degradation of the thermal contact in 
insulating powder samples in this temperature range. 

The previous spin-wave calculations^ reflect the gen- 
eral features of the measured spectrum, with a 10-20% 
accuracy for the observed gaps. The principal qualitative 
difference with our work is that our ESR measurements 
find two exactly degenerate magnon branches at k = 
in zero magnetic field, whereas Ref. [l4] predicts a finite 
splitting between them: A2 = 1.76 K and A3 = 1.93 K. 
In the next section we present the detailed theoretical 
calculations of the ESR spectra in Gd2Sn207, which not 
only yield the correct degeneracy of k = magnons but 
also show an overall improved agreement with the exper- 
imental data. 



est experimental temperature 0.45 K. 

The upper panel of Fig. [3] shows the resonance spectra 
for a sample with h mw _L H. These data have been briefly 
discussed before*^ The single resonance lines observed 
at various frequencies belong to two different branches, 
one of which is decreasing (line A) and the other one is 
increasing (line B) in field. The extrapolation of these 
lines to zero field gives the same gap value of 33.5 ± 
0.5 GHz (1.61 ± 0.02 K) for both branches, which points 
to an exact degeneracy of the corresponding magnetic 
excitations. No other resonance absorption was detected 
within the experimental accuracy of 0.5%. 

The absorption spectra obtained for a tilted direction 
of h mw are presented on the lower panel of Fig. [3] The 
intense sharp resonance line shown in Fig. [21 which arises 
in the ordered phase, appears to be a third branch with 
a gap of 85 ± 5 GHz (4.1 ± 0.2 K, line C). This mode 
decreases with increasing magnetic field and softens in 
the vicinity of the transition into the fully saturated 
phase. Extrapolating its field dependence to zero fre- 
quency, one can determine the value of the saturation 



IV. SPIN- WAVE THEORY 

The unit cell of the pyrochlore lattice contains four 
magnetic atoms. Their positions and the specific choice 
of the local axes adopted below is 

Pl = (0,0,0), ^ = -^(1,-1,0), 

p 2 =(0,ii), z 2 = -^(-1,-1,0), 

P3=(j,0,i), 33 = ^(1,1,0), 

p 4 =(I,i,o), z 4 = -i=(-l,l,0) (1) 

with yi = (0, 0, 1) and Xj = x Zj. The unit cell tetrahe- 
dra form a face-centered cubic structure with the transla- 
tional basis ax = (0, \, \), a 2 = (|,0, ±), a 3 = (±, ±, 0). 

In the following we will take into account all major 
magnetic interactions starting with the strongest nearest- 
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neighbor exchange: 



terms in the bosonic representation one obtains: 



H = J 



(2) 



where n,m denote unit cells and i,j = 1-4 indicate 
position inside cell. The equilibrium magnetic struc- 
ture determined for Gd2Sn2 07 in neutron diffraction 
experiments^ is a four-sublattice chiral spin-cross with 
Sj || Zj, which is also known as the Palmer-Chalker 
stated This state is one out of the many degenerate spin 
structures obeying the classical constraint: J2i &i = for 
every tetrahedron, which must be fulfilled by the ground 
state in the absence of single-ion anisotropy and dipolar 
interactions. 

Let us first calculate the excitation spectrum for 
the chiral spin-cross state neglecting all magnetic 
anisotropics. The local coordinate system is used for the 
spin operators Si = Sf 5q + Sf y, + Sf&i, where the three 
axes are defined above, Eq. (JT|). Since we are dealing with 
the k = magnetic structure, the local axes do not carry 
the cell index n. The Holstein-Primakoff transformation 
is used for the bosonic representation of spin operators. 
Since gadolinium spins S — 7/2 are large, we neglect in- 
teractions between magnons and always keep up to the 
quadratic terms in the bosonic Hamiltonian. With the 
chosen accuracy it is sufficient to write 



(3) 



Sf = S-4ai, Sf = \l-(al + ai ) , 



The obtained quadratic form is diagonalized by perform- 
ing consecutively Fourier and Bogoliubov transforma- 
tions. The excitation spectrum consists of four branches: 



£i, 2 (k) = 



_-:!. l(k) = 2 J.S'\/l-COS yCOSy 



(4) 



Two of them have zero energies everywhere in the Bril- 
louin zone, which reflects the infinite degeneracy of the 
nearest-neighbor exchange Hamiltonian on a pyrochlore 
lattice. Interestingly, the two other branches have purely 
two-dimensional dispersions in the harmonic approxi- 
mation. Various perturbations to the nearest-neighbor 
Heisenberg Hamiltonian © should produce a finite dis- 
persion of the two lowest branches, but the quasi 2D be- 
havior of e 3j4 (k) for the two higher-energy modes should 
survive in a certain range of parameters. 

As the next step we add the single-ion anisotropy to @ 
in the form of the lowest-order crystal-field harmonics: 



(5) 



rii being the four local anisotropy axes parallel to the 
principal cubic diagonals. Keeping again only quadratic 



(Si • n,) 5 



= s(a\ai + - 



(6) 



+ 



f(l ± 2V2i) +af(l±2y/2i 



where the upper (lower) sign corresponds to i = 1, 4 (2,3). 

Since we are interested in the ESR spectrum given by 
k = magnons, we simplify the following calculations by 
considering only uniform modes, avoiding hence a step 
with the Fourier transformation. The exchange Hamil- 
tonian projected onto the four-site magnetic unit cell is 
written as 

"Hi I J SN — 2{a\a\ + a\a 2 + a\a?, + a 4 a 4 ) 

+ a\(a 2 + a 3 ) + a\{ai + a 4 ) + a/ 3 (ai + a 4 ) 

+ a\{a 2 + a 3 ) - ai(a 2 + a 3 + 2a 4 ) 

— a 2 (2a 3 + a 4 ) — a 3 a 4 + h. c. , (7) 

where N is the number of unit cells in the lattice. The 
single-ion anisotropy is represented as 



H^/DaSN = a\ai+a\a 2 
1 



a' 3 a 3 



a 4 a 4 



(8) 



+ -(1 - 2V^i)(ffli + a\ + a 2 2 + a 3 2 ) + h. c. 

Diagonalization of quadratic forms with the help of a gen- 
eralized Bogoliubov transformation has been described 
many times in the literature and will not be repeated 
here. In the present case the diagonalization can be per- 
formed analytically yielding four magnon gap energies: 



Ai 



, A 2i3 = 2S^/2D a J/3 , A 4 = 2A 2 . (9) 



The lowest magnon branch remains gapless in spite of the 
single-ion term. This is in agreement with the analysis 
of Ref. [l(| which finds an infinite (but not extensive) de- 
generacy for the easy-plane pyrochlore antiferromagnet. 

The double degeneracy of the two intermediate 
magnon modes A 2 , 3 follows from the tetragonal sym- 
metry of the chiral spin-cross structure. Analysis of 
the eigenvectors of the Bogoliubov transformation iden- 
tifies corresponding oscillations with out-of-plane motion 
of only one pair of spins Si^S^ or S 2 ,S3 with opposite 
phases. The two modes transform into each other accord- 
ing to the 2D irreducible representation of the tetrago- 
nal point group. Extra interactions (dipolar, etc.) will 
not modify such a degeneracy as long as the chiral spin- 
cross structure remains stable. The two other modes Ai, 4 
correspond to predominantly in-plane motion of all four 
spins. 

For a simple numerical estimate we use the following 
parameters for Gd 2 Sn 2 7 : |0cw| = 2JS(S+1) « 8.6 K±I 
with J = 2/63 |0cw| = 0.27 K and D a sa 0.14 K. 17 
This yields the following values A 2 3 = 1.12 K and 
A 4 = 2.24 K, which are somewhat lower than the exper- 
imentally measured frequencies (1.61 and 4.1 K respec- 
tively). Introduction of further neighbor exchanges will 
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not modify the above results: (i) the third-neighbor ex- 
change J3 (the notation is taken from Ref. l3j) couples 
spins on the same sublattice and, consequently, does not 
contribute to the k = modes, (ii) the antiferromagnetic 
second-neighbor exchange J2 yields the same replacement 
J — ► J + 2 J2 in the expressions for Ocw and for the uni- 
form modes and does not, therefore, change the gaps. 

The next important perturbation to the exchange en- 
ergy (J2j) is the dipole-dipole interaction ; 14 i 18 



D J 

m,mj 



3(Sr, 



,)(S, 



,) 



with D = (gfiBS) 2 /a 3 , RJf TO being the vector linking 
two spins, measured in units of the lattice constant a, 
and r l J m — R^ m /|RJf m |. The strength of the dipolar 
coupling between two neighboring spins in Gd2Sn2C>7 is 
estimated as 



Et 



(gfi B S) 2 
{a/2V2f 



= WV2D = 0.605 K 



(11) 



where we have substituted a = 10.455A for the lattice 
constant^ The parameter E^ n is three-times smaller 



than the single-ion energy D a S 2 ~ 1.73 K, but the dipo- 
lar interactions still play an important role due to their 
long-range nature. From now on we define dimensionless 
parameters d a = D a /J = 0.516 and d = D/JS 2 = 0.008 
normalizing all interactions to the exchange constant J 
and the excitation energies to JS. 

Projecting (|10p onto the four-sublattice magnetic 
structure (S n i = Si) we obtain 



1 



-N 



E S?S?Dg 



(12) 



(10) 

where the dipolar matrix is given by 



<^E 



« 13 



da/3 ^( r nm) ( r nm) 



(13) 



In a cubic crystal the diagonal matrix elements are 
Dff ~ 8 a p and drop out from the equations on the 
equilibrium spin configuration and on the energies of the 
k = modes. 

The dipolar sums are straightforwardly evaluated using 
the Ewald's summation technique^ 
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lbn > — e 

G 



G 2 



G 2 /AQl iG-( Pi - Pj ) 



_2Qc_ 



-QoRl 



2Q c R? i R ii - 5 a 



E 

R 



erfc((3 c i? y ) 



c Q TDa DP " 

0af3 C,lt i3 n ij 



(14) 



where Q c ~ 1 is an arbitrary cutoff, G is a reciprocal lat- 
tice vector, erfc(a;) is the complementary error function, 
and Ry = R + p i — p 3 , R being an fee lattice vector. 
Summations in Eq. (114j) are performed over all R and G 
excluding G = 0. 

Cubic symmetry of the pyrochlore lattice leaves only 
three independent constants 
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(15) 



Evaluating numerically the corresponding expressions 
(fig)) we find ci = 17.92d = 0.143, c 2 = -34.09d = 
-0.273, and c 3 = -57.84d = -0.463. 

The dipolar energy (Tl2|) can be used to compare the 



relative stability of different q = magnetic structures. 
In particular, the dipolar contribution for the chiral spin- 
cross configuration (JXJ) , or the Palmer-Chalker state, is 
Ef/NS 2 = — 2ci + 2c3. The alternative (non-chiral) spin- 
cross structure, which is realized in Er2Ti207^i and can 
be obtained from the chiral spin cross with S2 — » — S2, 
S 3 — > —S3, has a higher dipolar energy E2/NS 2 = 2c\ — 
4c 2 + 2c 3 . 

To calculate the effect of the dipolar interaction on 
the magnon spectra we transform again to the local spin 
frame in Eq. (|12[) and bosonize spin operators using ([3]). 
The dipolar matrix elements in the rotating coordinate 
system "D"? are expressed via the laboratory frame ma- 
trix D«f by 

V°f = e^e» D%, (16) 

where e.i a are the local basis vectors, see Eq. (fT]). 
The obtained quadratic form of bosonic operators is 
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''•dip 



E 



XX 



+ \{V%-Vyy){a i a 3 +a\cJ j )- i ~{V 



xy 
ij 



(17) 



The summation is performed over all sublattice pairs. 
Skipping the straightforward algebra behind the substi- 
tution (|16p and the subsequent Bogoliubov transforma- 
tion, we present the final results for the magnon energies. 
The two degenerate modes have the energy: 

£2,3 = \ [(2do - 3c 3 ) (4 + ci + ca) - 4c 3 d a ] , (18) 

while the energies of the two other modes are given by 
(positive) roots of 

3£ 4 -2£ 2 [6(ci-c 2 )(4 + ci+c 2 ) (19) 
+ (2d a -3c 3 )(8 + 5 Cl -c 2 -4c 3 )] 
+ 8(ci - ca) [2d a - 3c 3 + 2(ci - ca)] 
x [(4 + ci + c 2 )(2d - 3c 3 ) - 4c 3 d Q ] = . 

Restoring the scaling parameter JS and using the follow- 
ing values for the microscopic parameters: J = 0.27 K, 
d a = 0.516, and d = 0.008 we obtain for the gaps: 
Ax = 1.24 K, A 2i3 = 1.77 K, A 4 = 4.51 K, which are 
already quite close to the experimentally measured val- 
ues. An even better correspondence between the theo- 
retical and the experimental results for the three upper 
branches (shown by closed squares on the H = axis of 
Fig. is achieved for a slightly modified set of micro- 
scopic constants: J = 0.25 K, D a = 0.13 K (d a = 0.52), 
and E d = 0.55 K (d = 0.008). The lowest gap in this 
case is equal to Ai = 1.13 K, which is very close to the 
estimate obtained by fitting the specific heat curve. 



V. DISCUSSION 

In the preceding section we have demonstrated a good 
agreement between the measured ESR spectrum and the 
spin-wave calculations performed in zero magnetic field. 
The two degenerate resonance branches correspond to a 
peculiar type of spin motion: oscillation of a spin plane 
with respect to two orthogonal in-plane axes. Such modes 
are excited by the perpendicular component of the mi- 
crowave field h mw 1 H as is indeed observed for the 
lines A and B. 

Generalization to weak magnetic fields is straightfor- 
wardly done with the "hydrodynamic" approach, 22 which 
is valid once the exchange interactions are significantly 
stronger than magnetic anisotropics and field. Simple 
calculation analogous to those performed in Ref. [24j yield 
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FIG. 4: (Color online.) Frequency-field diagram of the reso- 
nance spectrum observed on a powder sample of Gd2Sn207 
at T — 0.45 K; O and © correspond to the spectra obtained 
at h mw _!_ H and h mw ZH respectively; ■ mark the gap val- 
ues calculated by the spin- wave theory (SWT); the solid lines 
represent the solution of Eq. (|20p for a doublet mode with 
H || [001], the dashed lines are linear fits to the two compo- 
nents of spectral line E with g — 2, and the dash-dotted line 
corresponds to a g = 2 paramagnet. The transition to a spin- 
polarized phase at H s = 53.0 kOe is marked by a vertical line. 
The unreachable low-frequency range < 25 GHz is shaded 
in grey. 

the following cubic equation for eigen-frequencies: 

(20) 

- 7 V [v 2 H 2 ~ v\B\ - vlH 2 ^j = 0, 

where V\ and ^ 2 are resonance frequencies in zero field, 
H\\ and H± are magnetic field components with respect 
to the tetragonal axis, 7 = g/is/27rfi is the electronic 
gyromagnetic ratio (g = 2.0). In gadolinium stan- 
nate the magnetic anisotropies are comparable to the 
nearest-neighbor exchange interaction. Moreover, the ex- 
change structure is soft, i.e., infinitely degenerate, and 
the anisotropies play a decisive role in stabilizing the ob- 
served magnetic structure. This restricts applicability 
of the hydrodynamic theory, which is used here only to 
indicate a plausible behavior. 

The field evolution of the spectrum is summarized on 
the frequency-field diagram presented in Fig. 3J The two 
degenerate modes appear to be split by the magnetic 
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field into decreasing (Line A) and increasing (Line B) 
branches. The points in Fig. 0] mark the maximums of 
the ESR absorption which, for a powder sample, cor- 
respond to one of the outermost field orientations with 
respect to the crystal axes. The frequency-field depen- 
dence of these maximums for the two degenerate modes 
is satisfactorily fitted by formula (f^Uj) for H | [001] if one 
takes hv2 = A2.3 (solid lines on Fig. Q|, while the third 
calculated branch is field independent and set to zero. 
When the orientation of the external field is changed 
from H || [001] -> H _L [001], both branches shift to 
higher fields resulting in the overextended right wings of 
the absorption lines observed in the experiment. Hence, 
in spite of the "softness" of the exchange structure in a 
Heisenberg pyrochlore magnet, the spin plane oscillations 
are not strongly affected by quasi-local modes. 

The third spectral mode C is excited only by a paral- 
lel microwave field component /i mw || H , which indicates 
that it is not a uniform oscillation of the spin plane, but 
rather an antiphase motion of spins of the cross (in-plane 
or out-of plane). This mode should soften at the satu- 
ration point H s (as observed in the experiment) and, on 
further increase of H, should develop a gap, which prob- 
ably corresponds to the increasing line D of the spectrum 
above H s . The other two resonance branches (the dou- 
blet line E) observed at H > H s have linear field depen- 
dences distinctive for quasi-local soft modes in the spin- 
polarized phase.- 3 They were also observed in Gd2Ti2C>7 
but unlike that case, one of them has almost zero en- 
ergy at H = H s , which corresponds to the softening 
of one of the excitation branches at the antiferromag- 
netic wave- vector q = near the second-order transition. 
(The ordering wave- vector in Gd2Ti2C>7 is different from 
q = and, therefore, all ESR modes remain finite at 
H = H s .) The spin- wave calculation of the high-field 
magnon spectrum and its comparison with the observed 
results can provide a useful information on the parame- 
ters of the the spin-Hamiltonian for the two pyrochlore 
materials. Also, an interesting question remains about 
possible phase transitions in Gd2Sn2C>7 at intermediate 
fields < H < H 8 . Such transitions, if any, should be 
quite sensitive to the field orientation, which hinders their 
observation on the available powder samples. 

The measured resonance modes and their perfect 
agreement with the spin-wave calculations point at the 
conventional magnetic ordering in Gd2Sn2C>7 below 1 K. 
Note, that an explanation of the persistent spin dynam- 
ics observed by the local probesi^ remains an open is- 
sue: gapped magnons with energies above 1 K cannot 
contribute to the muon spin relaxation in the 50 mK 
range. Still, properties of the sister material Gd2Ti2 07 
are strikingly different from the more conventional stan- 
nate. This concerns not only the difference in the low- 
temperature asymptotes for the specific heat^ but also 
their ESR spectra as illustrated in Fig. [5] A typical 
absorption curve for Gd2Sn207 consists of several spec- 
tral lines with almost zero background. In contrast, the 
recorded absorption curve for Gd2Ti207 demonstrates a 
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FIG. 5: (Color online.) Absorption spectra recorded at nearly 
the same frequencies in Gd2Sn207 powder sample (squares) 
and in Gd2Ti2C>7 single crystal at H || [111] (circles); the ar- 
row marks the transition in Gd2Ti2C>7 at H = H c \ ~ 30 kOe. 



broad intense non-resonant anomaly with a maximum in 
zero field. This additional absorption develops simulta- 
neously with the main spectral lines below Tjv, rapidly 
decreases with increasing magnetic field and fully dis- 
appears stepwise at H = H c \ ~ 30 kOe marked by an 
arrow in Fig. [5] The above findings provide an evidence 
for additional "magnetic degrees of freedom" in the ti- 
tanate, which exist down to low temperatures and are 
suppressed by an external magnetic field. 

Additional low-energy excitations in Gd2Ti207 may 
result from the complexity of magnetic ordering in this 
material. The neutron diffraction experiments 2 suggest 
a multi-A: spin structure in the titanate, which may 
also lead to multiple magnetic domains related to dif- 
ferent combinations ofq = (1/2, 1/2, 1/2) and equivalent 
wave-vectors. Excitations in the domain walls and their 
pinning by crystal defects can produce additional low- 
temperature spin dynamics. Such a residual dynamics 
(i) should be absent in the q = ground-state of the 
stannate and (ii) exhibit a significant sample dependence 
as was indeed found from the comparison of /iSR data 
obtained on single crystal^ and powder samples^ The 
low-field domain structure can be further eliminated by 
a magnetic field, which, for example, selects at H > H c \ 
a unique ordering wave- vector. An extra argument in fa- 
vor of such a scenario is that the ESR spectrum of the 
titanate is significantly transformed at H = H c \ becom- 
ing similar to the spectrum of the stannate^ Namely, 
the line C of the spectrum, which is traced in the stan- 
nate in the whole field range < H < H s , appears in 
the titanate only at H > H c i . High-field neutron diffrac- 
tion measurements in Gd2Ti207 together with evolution 
of the low-T asymptote in the specific heat under mag- 
netic field should provide a valuable check for the above 
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scenario. 

In summary, the study of the magnetic resonance prop- 
erties of the pyrochlore gadolinium stannate reveals three 
gapped resonance modes in the ordered phase, two of 
them being exactly degenerate at zero external magnetic 
field. The spin-wave theory, which takes into account 
the nearest-neighbor exchange, the single-ion anisotropy 
and the dipolar interactions, demonstrates a very good 
agreement with the experiment using known values of 
the microscopic magnetic parameters. The lowest gap 
value predicted by theory lies beyond the experimental 
frequency range and cannot be directly observed. Nev- 
ertheless, it roughly agrees with an estimate made from 
the exponential decrease of the low temperature specific 



heat. 
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